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ABSTRACT 

A  packet  switching  network  has  the  desirable  feature  of 
rapidly  handling  short  (bursty)  messages  of  the  type  often 
found  in  computer  communication  systems.  In  evaluating  packet 
switching  networks,  the  average  time  delay  per  packet  is  one 
of  the  most  important  measures  of  performance. 

The  problem  of  message  routing  to  minimize  time  delay  is 
analyzed  here  using  two  approaches,  called  "successive  satura¬ 
tion"  and  "max-slack",  for  various  traffic  requirement  matri- 
cies  and  networks  with  fixed  topology  and  link  capacities. 
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INTRODUCTION 


A.  THE  PACKET  SWITCHED  NETWORK  CONCEPT 

A  new  technique  for  data  communications  that  has  evolved 
over  10  years  is  called  PACKET  SWITCHING.  In  general,  commu 
nication  networks  may  be  conveniently  divided  into  three 
types:  CIRCUIT  SWITCHING,  MESSAGE  SWITCHING  and  PACKET 

SWITCHING.  3oth  message  and  packet  switching  uses  a  tech¬ 
nique  known  as  store  and  forward  transmission. 

A  circuit  switching  network  provides  service  by  setting 
up  a  total  path  of  connected  lines  from  the  source  to  the 
destination  of  the  call.  This  complete  circuit  is  set  up  by 
a  special  signaling  message  that  threads  its  way  through  the 
network,  seizing  channels  in  the  path  as  it  proceeds.  After 
the  path  is  established,  a  return  signal  informs  the  source 
that  daca  transmission  may  proceed,  and  all  channels  in  the 
path  are  then  used  simultaneously. 

The  entire  path  remains  allocated  to  the  transmission, 
and  only  when  the  source  release  the  circuit  will  all  these 
channels  be  returned  to  the  available  pool  for  use  in  other 
paths.  Circuit  switching  is  the  common  method  for  telephone 
systems. 

In  message  switching,  only  one  channel  is  used  at  a  time 
for  a  given  transmission.  The  message  first  travels  from 
its  source  node  to  the  next  node  in  its  path,  and  when  the 
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entire  message  is  received  at  this  node,  then  the  next  step 
in  its  journey  is  selected-  If  this  selected  channel  (link; 
is  busy,  the  message  waits  in  a  queue,  and  finally  when  the 
channel  becomes  free,  transmission  begins.  Thus,  the  message 
"hops"  from  node  to  node  through  the  network  using  only  one 
channel  at  a  time,  possibly  queueing  at  busy  channels,  as  it 
is  successively  stored  and  forwarded  through  the  network. 

Packet  switching  is  basically  the  same  as  message  switch¬ 
ing  except  that  the  messages  are  decomposed  into  small  equal 
pieces  called  packets,  each  of  which  has  a  minimum  length. 

These  packets  are  numbered  and  addressed  and  make  their  way 
through  the  net  independently  of  each  other.  Thus,  many 
packets  of  the  same  message  may  be  in  tranmission  simulta¬ 
neously,  giving  one  of  the  main  advantages  of  packet  switching. 

With  packet  switching  systems,  information  is  exchanged 
in  the  form  of  short  packets.  A  packet-switched  r.  2twork  can 
handle  several  different  types  of  traffic  concurrently.  These 
include  HIGH-THROUGHPCJT  traffic,  for  example,  the  transmission 
of  large  data  files  between  computers,  for  which  accuracy  and 
high  average  data  speed  are  the  most  important  performance 
requirements;  LOW-DELAY  traffic,  for  example,  interactive 
communication  between  a  person  at  a  terminal  and  a  remote 
computer,  for  which  accuracy  and  low  average  message  delay  are 
important;  and  REAL-TIME  traffic,  for  example,  packetized 
speech  for  which  the  performance  of  circuit-switched  connec¬ 
tion  must  be  approached  by  maintaining  a  relatively  constant 
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data  speed,  but  for  which  extreme  accuracy  is  not  important 
owing  to  the  redundancy  of  the  information. 

The  packet  switched  network  is  designed  primarily  for 
computer  to  computer  communication.  It  has  a  much  more 
rapid  response  which  matches  the  internal  behavior  of  compu¬ 
ters  and  handles  information  in  much  the  same  way  as  does  a 
computer.  At  the  same  time  it  can  readily  match  the  speed 
of  attached  computers  to  that  of  the  terminal  users,  by  virtue 
of  its  internal  storage. 

The  prime  purpose  of  store  and  forward  packet  switching 
is  to  enable  communications  resources  to  be  used  effectively 
and  in  such  a  way  that  they  may  be  shared  by  many  users  opera¬ 
ting  in  an  intermittent  fashion,  giving  each  user  a  rapid 
response  from  the  communication  network  just  at  the  instant 
when  this  is  required. 

If  there  is  need  for  transmitting  a  long  continuous  stream 
of  data,  then  a  circuit  switched  connection  makes  good  sense. 
On  the  other  hand,  if  the  data  flow  is  bursty,  then  some  form 
of  resource  sharing  can  be  used  to  great  advantage;  packet 
switching  is  an  effective  choice  here. 

Since  packets  are  stored  as  they  pass  through  switching 
nodes,  it  is  possible  to  conduct  speed,  format  and  code  con¬ 
version  during  the  switching  process.  Another  feature  of 
packet  switching  is  its  ability  to  adaptively  select  good 
paths  for  packet  journeys  as  a  function  of  the  network 
congestion. 
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3esides  providing  small  network  delays,  packet  switching 
has  the  desirable  feature  of  rapidly  handling  small  messages 
in  spite  of  the  presence  of  long  messages  that  may  be  in 
transport  at  the  same  time,  this  is  because  of  the  decomposi¬ 
tion  of  long  messages  into  packets.  Another  useful  property 
of  this  decomposition  is  that  the  nodal  storage  requirement 
is  reduced. 

In  evaluating  packet  switching  networks,  the  delay, 
throughput,  cost  and  reliability  are  important  measures. 

Theoretical  studies  have  been  directed  to  the  queueing  and 
network  flow  problems  in  general  and  more  specifically  to  such 
problems  as  delay  analysis,  route  assignment,  topological 
design  and  flow  control,  etc. 

The  topic  chosen  in  this  study  is  centered  around  optimal 
flow  and  minimum  capacity  assignment  and  delay  analysis  for 
packet-switched  networks.  The  networks  under  consideration 
have  a  relatively  complicated  structure  with  a  large  number 
of  source-destination  node  pairs. 

B.  THE  ROUTING  PROBLEM 

A  message  routing  procedure  is  merely  a  decision  rule 
that  determines  the  node  a  message  will  next  visit  in  its 
path  through  the  network.  The  objective  of  the  routing 
procedures  is  to  transport  packets  on  a  minimum  delay  path 
from  source  to  destination. 

In  discussing  routing  policies  for  networks,  an  important 
distinction  must  be  made  between  static  and  adaptive  policies. 
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This  distinction  depends  on  the  environment  in  which  a  policy 
is  designed  to  operate.  If  network  topology  is  not  subject 
to  changes  (due  to  failure,  modifications,  growth)  and  traffic 
inputs  are  stationary,  then  the  optimal  routing  solution  is  a 
static  solution  consisting  of  a  set  of  fixed  paths  between  all 
node  pairs.  The  traffic  between  each  source  and  destination 
pair  may  be  distributed  on  several  paths  simultaneously  in 
well  defined  proportions,  where  the  proportions  are  fixed  in 
time . 

In  a  real  network  environment,  however,  the  topology 
changes  with  time  and  user  traffic  requirements  tend  to  fluc¬ 
tuate  more  or  less  rapidly.  To  minimize  delay  it  is  then 
necessary  to  implement  an  adaptive  routing  policy  that  can 
react  and  adjust  to  various  changes.  The  adaptivity  of  a 
policy  can  be  measured  in  terms  of  its  response  time.  A 
reasonable  procedure  is  to  use  a  periodically  refreshed 
static  routing  solution.  Indeed,  there  is  a  continuum  of 
solutions  between  the  two  extremes,  static  and  adaptive, 
characterized  by  different  response  times  and  used  for 
different  applications. 

Beside  their  use  in  operational  networks,  static  routing 
policies  find  an  important  application  in  the  network  design 
process.  During  this  process,  for  a  given  traffic  pattern  a 
prediction  of  the  throughput  and  delay  performance  of  a 
given  topolocy  is  needed.  The  routing  policy  clearly  has  a 
major  impact  on  such  performance.  Most  routing  implementa¬ 
tions  are  adaptive,  and  unfortunately  the  analysis  of  adaptive 
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routing  policies  is  an  extremely  difficult  task.  To  simplify 
the  problem,  the  traffic  pattern  is  usually  approximated  with 
a  stationary  pattern,  and  the  routing  policy  with  a  static 
routing  policy. 

1 .  Static  Routing 

a.  Routing  Table  Representation 

A  static  routing  policy  is  represented  by  a  set 
of  routing  tables,  one  for  each  node,  indicating  how  the 
traffic  arriving  at  that  node  must  be  routed  on  the  outgoing 
links,  depending  on  its  final  destination.  The  routing  table 
for  node  i  is  an  Nx  L  matrix  Pi(n,k) ,  where  N  is  the  number 
of  nodes  in  the  network  and  Pi(n,k)  is  the  fraction  of  traffic 
directed  to  node  k  which,  upon  arrival  at  node  i,  is  routed 
through  neighbor  node  n  (Fig.  1.1).  By  the  definition. 

Pi  (j,k)  =  1 

Where  n  is  the  number  of  neighbors  of  node  i. 

The  actual  distribution  of  the  incoming  traffic 
to  outgoing  links  may  be  done  randomly  using  as  probability 
weights  the  values  Pi(j,k).  If  for  any  node  i  there  is  only 
one  permissible  outgoing  link  from  i  to  any  destination  k,  then 
Pi(j,k)  =  1  for  neighbor  node  j. 

0  otherwise 

Such  a  routing  policy  is  referred  to  as  a  single  path 
routing  policy,  since  only  one  path  is  used  from  any  node  i 
to  any  destination  node  k.  In  general,  the  optimal  static 
routing  solution  is  a  multipath  solution,  allowing  for  the 
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simultaneous  use  of  several  routes  in  order  to  minimize 
delays.  The  routing  table  representation  may  also  be  used 
for  adaptive  policies. 


F i  (-  f  k J  j 

p  in , 

_ 

Figure  1.1.  Routing  Table. 


b.  Problem  Formulation 

In  packet  switched  network,  messages  are  seg¬ 
mented  into  small  packets,  and  each  packet  travels  from  its 
source  to  its  destination  via  a  set  of  intermediate  nodes. 
While  awaiting  transmission  to  the  next  node  in  the  route, 
the  packet  must  be  stored  until  the  link  to  that  node  is 
free.  Thus,  at  each  node  there  are  several  queues,  one  for 
each  output  channel. 
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Packet  flow  requirements  between  nodes  arise  at 


random  times;  therefore,  link  flows,  queue  length,  and  packet 
delay  are  random  variables.  A  static  routing  problem  can  be 
defined  as  the  problem  of  finding  the  static  routing  policy 
which  optimizes  the  flow  and  minimizes  the  average  time  delay. 

Considering  the  relationship  between  routing  policy 
and  link  flows,  the  problem  can  be  formulated  for  solution  of 
optimal  flow  and  minimum  capacity  with  respect  to  the  given 
requirement  as  follows, 
given:  (a)  topology 

(b)  requirement  matrix  R 
minimize:  alpha 

subject  to:  (a)  conservation  of  flow 

(b)  flow  <_  capacity  x  alpha 
Where  alpha  is  a  nonnegative  scale  factor. 

The  mathematical  model  and  solution  techniques  are 
discussed  in  detail  later. 

2 .  Adaptive  Routing 

The  adaptive  routing  problem  consists  of  defining  a 
procedure  that  dynamically  updates  routing  tables  according 
to  changes  in  the  network.  This  procedure  will  include  the 
acquisition  of  the  status  of  all  neighboring  nodes,  which  a 
packet  associated  with  the  source  and  destination  node  pair 
could  visit  next  after  being  processed  at  the  current  node. 

The  neighboring  nodes  are  those  with  direct  channels  linking 
to  a  given  node. 


The  structure  of  the  routing  table  and  the  ways  the 
information  contained  therein  are  to  be  used  vary  according 
to  the  individual  routing  procedure.  For  example,  if  single 
path  routing  is  used,  only  one  path  links  each  source  and 
destination  node  pair;  therefore,  only  a  single  neighboring 
node  associated  with  each  node  pair  is  listed  in  the  routing 
table . 

The  design  of  routing  is  a  process  of  determining  the 
structure  of  routing  tables  and  specifying  the  procedures  of 
using  them.  The  choice  of  routing  depends  on  many  factors 
such  as  network  topology,  available  facilities,  throughput 
and  delay  requirements,  and  so  on.  Adaptive  routing  tech¬ 
niques  are  ruled  out  in  this  study. 
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II .  OPTIMAL  FLOW  AMD  MINIMUM  CAPACITY  ASSIGNMENT 

The  optimal  flow  is  the  best  selection  of  paths  between 
source  and  destination,  given  the  traffic  requirements  and 
the  network  configuration.  The  definition  of  best  paths  may 
vary  depending  on  the  nature  of  the  traffic.  In  this  study, 
the  best  paths  will  be  regarded  as  the  paths  of  minimum 
"distance",  where  distance  can  be  interpreted  in  various  ways, 
for  instance  as  time  delay  along  the  links. 

To  investigate  average  time  delay  of  the  networks,  two 
different  approaches  are  discussed  using  classical  optimiza¬ 
tion  procedures . 

A.  NETWORK  MODELING 

In  a  distributed  packet-switched  network,  packets  are 
transmitted  over  a  collection  of  switching  nodes  in  tandem. 

In  modeling  a  packet  switched  network,  the  objective  is 
usually  to  optimize  the  total  flow  of  the  networks,  to 
minimize  the  capacity  of  each  link,  or  to  maximize  the 
utilization  of  the  network.  The  problem  may  be  formulated 
as  follows. 

given  topology 

requirement  matrix  R 

capacity  of  each  link  (arbitrary) 

minimize:  alpha 

subject  to:  (1)  conservation  constraints. 
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>  Fc (k)  -  lFi(k)  =  Ri (k) 

outgoing  incoming  input 

links  links  at  node  i 

from  node  i  to  node  i  for  node  k 

(2)  capacity  constraints 
F-  (k)  £  alpha  x  CL'. 

(3)  positive  constraints 
I  Ft  (k)  _>  0 

Where  Ft (k)  is  the  flow  on  link  destined  for  node  k.  We 
next  look  at  each  constraint  further  in  detail. 

1 .  Conservation  Constraints 

At  each  node  i,  the  total  outgoing  packet  flows 
destined  for  node  k  equal  the  incoming  flows  plus  the  input 
to  be  transmitted  to  node  k. 

I  Fi  (k)  -  If  UK)  =  Ri(k) 

outgoing  incoming  input 

links  links  at  node  i 

from  node  i  to  node  i  for  node  k 

whe  re 

k  =  destination  node  (k  =  1,2,3...N). 

N  =  number  of  nodes. 

L  =  number  of  links. 

i  =  source  node  (i  =  k,  i  =  1,2,3. . .N) . 
l  =  link  between  two  nodes  (i=  1,2,3...L). 

Fi(k)  =  flow  on  link  destined  for  node  k. 

Ri(k)  =  input  destined  for  node  k  at  node  i. 
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2 .  Capacity  Constraints 

The  total  flow  destined  for  node  k  on  a  certain  link 
should  be  less  than  or  equal  to  the  capacity  of  the  link  l. 
Since  linear  programming  algorithms  have  difficulties  dealing 
with  inequalities,  slack  variables  are  added  to  the  constrains 
to  absorb  unused  resources  and  thus  to  force  equalities  to 
appear. 

When  slack  variables  are  introduced,  the  problem  in 
standard  form  becomes. 

Fc(k)  -  alpha  X  Cc  +  Sc  =  0 

k=l 

where 

Ct  =  capacity  of  the  link. 

alpha  =  constant  to  optimally  minimize  the  link 
capacity  and  total  flow  on  link. 

S i  -  slack  variables. 

3 .  Positive  Constraints 

If  we  consider  the  network  to  have  unidirectional 
links,  the  bidirectional  flows  between  two  nodes  are  divided 
into  two  one-way  links  in  opposite  direction,  and  the  flow  on 
each  link  is  positive. 

Ft  (k)  _>  0 

Using  matrix  notation,  the  model  simply  may  be 
expressed  as  a  linear  programming  problem. 

minimize:  alpha 

subject  to:  AX  =  R 

x  >0 
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B.  MODEL  PROGRAMMING 

In  order  to  use  a  readily  available  Mathematical  Program¬ 
ming  System  such  as  MPSIII,  which  can  solve  linear  programming 
problems  with  up  to  4000  rows  and  theoretically  with  an  un¬ 
limited  number  of  variables,  the  major  programming  problem 
with  large  scale  networks  is  how  to  generate  the  input  data 
for  the  MPSIII  program.  In  other  words,  how  do  we  generate 
the  matrix  form  of  the  model. 

1.  Node  to  Link  Incidence  Matrix 

For  any  network  processing  to  take  place  in  the  compu¬ 
ter,  the  network  structure  must  be  represented  in  some  machine 
understandable  form.  A  wide  spectrum  of  different  represen¬ 
tation  methods  are  available.  Most  of  these  are  based  on  an 
incidence  matrix  scheme.  A  matrix  form  is  convenient  for  the 
analysis  of  networks  since  the  resulting  matrices  are  in  a 
format  suitable  for  mathematical  analysis. 

The  network  of  Fig.  2.1  will  be  used  to  illustrate  the 
method.  The  node  to  link  incidence  matrix  E(i,c)  describes 
the  links  connected  to  the  node  such  that 

E(i,i)  =  '”  I  if  link  c  is  outgoing  from  node  i. 

-1  if  link  i  is  incoming  to  node  i. 

0  if  link  t  is  not  connected  to  node  i. 


The  incidence  matrix  for  input  at  node  i  destined  for  node  k 
is  a  modification  of  the  node  to  link  incidence  matrix  E(i,£) 
according  to  the  following  rules: 


E(i,c)  is  a  3x5  matrix  and  Ek(i ,1)  is  a  2x5  matrix  in  size. 


Figure  2.1.  Illustration  Network. 

2 .  Matrix  Representation 

Matrix  representation  is  a  convenient  tool  for 
programming.  The  complete  matrix  representation  for  the 
illustration  network  is  shown  in  Appendix  A. 

The  general  model  matrix  AX  =  R  may  be  represented 
as  follows: 


A 

= 

R 

El 

Ft  (1) 

Ri  (1) 

E2 

Ft  (2) 

Ri  (2) 

S3 

Ft  (3) 

. 

-(3) 

: 

Ek 

' 

• 

Ft  (k) 

• 

Ri  ( k ) 

11  12  13  .  .  .  Ik  -  Ql 

Ft  (k) 

°' 

L 

st 

Alpha 

0 
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where 


Ek  =  (N-l)  x  L  incident  matrix. 

Ik  =  L  x  L  identity  matrix  with  zero  elements 

corresponding  to  zero  columns  of  the  Ek  matrix. 

I  =  L  x  L  identity  matrix. 

Ri(k)  =  requirement  matrix  which  is  the  inputs  to  be 

sent  from  source  node  i  to  destination  node  k. 

The  size  of  the  model  in  matrix  representation  13; 

A  =  ((N-l)  x  N  +  L)  x  ((N-l)  x  L  +  1) 

X  =  ((51+1)  x  L  +  1)  x  1 

R  =  ((N-l)  x  N  +  L)  x  1 

For  example,  the  network  consisting  of  20  nodes  and 
40  links  has  the  matrix  which  is  A  =  420  x  841,  X  =  341  x  1 
and  R  =  420  x  1  in  size. 

3 .  Programs 

To  generate  the  MPSIII  input  data  and  to  solve  the 
model,  requires  two  programs  which  are  listed  at  the  end  of 
this  study  (Appendix  C  and  D) . 

The  first  one  is  the  data  (model)  generation  program 
(Appendix  C)  written  in  the  Fortran  language.  The  output  of 
this  program  uses  special  notation  to  designate  the  flow  of 
the  link  connecting  node  pairs,  destination  node  of  the  flow, 
flow  variables,  slack  variables  and  input  at  the  node.  These 
notations  are  composed  of  a  letter  which  is  one  of  C,  L,  X  or 
S  and  7  numerical  digits. 


/ 
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The  notation  Cl,  Ll ,  Xc  ,  and  Sl  imply  conservation 
constraints,  capacity  constraints,  flow  variables,  and  slack 
variables  respectively.  The  rest  of  the  six  numbers  followin 
Cl,  Ll,  Xl  and  St  indicates  link  connecting  node  pair  with 
first  4  digits  and  destination  node  or  slack  variable  number 
with  last  2  digits.  For  example,  the  notation  C1OOO201  indi¬ 
cates  the  conservation  constraint  row  which  is  the  input  at 
node  2  destined  for  node  1,  L1Q2Q100  the  link  capacity  row  fo 
the  link  connecting  between  node  2  and  node  1,  X1O201O1  the 
flow  variable  on  the  link  connecting  between  node  2  and  node 
1  destined  for  node  1  and  S1000001  the  first  slack  variable. 

The  second  program  is  the  MPSIII  program  developed  by 
Management  Science  Systems,  Inc.  It  provides  elaborate  con¬ 
trol  language  for  the  formulation  of  solution  strategies  for 
mathematical  programming  problems. 

Individual  instructions  in  this  language  bring  in 
quite  elaborate  sections  of  code  designed  to  execute  the  step 
as  efficiently  as  possible. 

Using  the  special  notation  discussed  above,  there  is 
a  limitation  on  the  size  of  the  problem  which  can  be  solved, 
because  data  names  must  consist  of  8  characters  maximum. 

This  means  that  the  maximum  network  size  must  be  less  than 
100  nodes,  if  the  computer  capacity  allows. 

C.  SOLUTION  TECHNIQUES 

As  mentioned  before  ,  the  classicial  optimization  proce¬ 
dures  to  solve  the  investigation  model  provide  solutions 
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which  are  the  optimal  paths  and  minimum  link  capacities  for 


the  network.  Different  optimal  solutions  may  be  expected 
depending  on  different  methods  of  optimization.  The  two 
different  approaches,  each  involving  iteration  of  the  linear 
programming  solution  procedure  are  discussed  below,  and  illus¬ 
trated  in  terms  of  one  example  network  problem.  These 
approaches  are  based  on  the  characteristics  of  linear  pro¬ 
gramming  problems. 

1 .  Successive  Saturation  Approach 

The  standard  procedure  of  the  MPSIII  program  listed 
at  the  end  of  this  study  (Appendix  D  and  E)  produces  an 
"optimal"  solution  for  alpha  from  the  output  data  of  the  data 
generation  program  in  the  first  run.  But  the  corresponding 
flow  solution  is  not  a  good  solution  in  terms  of  the  average 
time  delay  in  the  packet-switched  network,  when  it  is  ana¬ 
lyzed  by  the  Kleinrock's  delay  analysis  model  (Ref.  4) .  The 
problem  is  that  only  the  flow  thru  the  saturated  links  is 
optimum. 

In  order  to  reduce  the  time  delay.  Iteration  runs  are 
attempted  until  all  links  are  saturated.  The  detailed  itera¬ 
tion  method  for  the  example  network  (Fig.  3.1)  which  has  5 
nodes  and  6  links,  each  with  original  link  capacities  10  (to 
simplify  the  problem)  is  discussed  below. 
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Figure  3.1.  Example  Network. 

The  ouput  of  MPSIII  program  is  divided  into  row 
section  and  column  section  (Ref.  9  and  Appendix  D  and  E) . 

The  objective  function  value  (which  is  the  same  as  the  alpha 
value)  and  the  links  saturated  with  the  alpha  value  in  the 
capacity  constraint  rows  indicated  with  the  first  letter- 
notation  L  can  be  read,  and  the  flows  on  each  link  and  the 
alpha  value  in  column  section.  In  the  example  network,  for 
the  first  iteration  the  saturation  level  (alpha  value)  is 
0.25,  and  the  saturated  links  are  L1010300  (link  3)  and 
L1040300  (link  7).  Their  activity  and  slack  activity  levels 
are  all  zero,  but  their  dual  activity  level  is  not  (pp  48)  . 

For  the  next  iteration,  the  changed  input  datas  are 
that  the  original  capacity  multiplied  by  the  saturation  level 
of  first  run  are  moved  from  the  alpha  column  to  the  right 
hand  side  (RHS)  column.  In  the  example,  the  new  link  capaci¬ 
ties  2.5  (0.25  x  10)  are  moved  to  RHS  (pp.  50  ),  and  the 
problem  is  solved  again.  Repeating  this  iteration  procedures 
until  all  links  are  saturated  and  moved  to  RHS,  the  objective 
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function  value  becomes  zero.  When  the  objective  function 
value  is  zero,  the  final  flow  solution  is  the  desired  solution 
of  this  successive  saturation  approach.  In  the  example  case, 
the  six  iterations  are  required  to  solve  the  problem.  These 
whole  procedures  except  programs  and  input  datas  from  2nd  run 
for  the  example  network  are  listed  at  the  end  of  this  study 
(Appendix  D) . 

2.  Max- Slack  ADoroach 

To  simplify  the  solution  of  the  above  successive  satu¬ 
ration  approach,  another  approach  called  the  max-slack  approach 
is  considered.  It  requires  only  two  iteration  of  the  program, 
but  is  otherwise  quite  similar.  Like  in  the  successive  satu¬ 
ration  approach,  after  the  first  run,  all  saturated  link 
capacities  are  multiplied  by  the  alpha  value  of  the  first  run 
are  moved  to  the  R HS,  and  in  addition,  all  slack  variable 
columns  are  summed  as  the  objective  function  in  the  data  card 
for  the  second  run. 

In  the  example  network,  eight  slack  column  variables 
are  added  on  the  objective  rows,  and  all  new  capacities  2.5 
(0.25  x  10)  are  moved  to  RHS .  This  is  clearly  shown  on  the 
output  of  control  language  "picture”,  and  the  program,  input 
data  and  the  output  are  also  listed  at  the  end  (Appendix  E) . 

D.  DELAY  ANALYSIS 

The  average  time  delay  of  a  network  is  the  average  time 
a  packet  spends  in  the  network  traveling  from  its  source  to 
its  destination  nodes.  To  obtain  a  tractable  expression  for 
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average  time  delay  T  in  terms  of  the  capacity  of  link  it 
is  necessary  to  make  simplifying  assumptions.  Packet  lengths 
are  assumed  to  be  exponentially  distributed,  and  are  re-chosen 
with  statistical  independence  at  each  node  to  form  a  Poisson 
process . 

These  assumptions  known  as  the  'Independence  assumption', 
introduced  by  D.  Kleinrock  (Ref.  4) ,  make  each  queuing  prob¬ 
lem  independent.  In  the  analysis,  each  queue  in  a  packet 
network  is  assumed  to  be  an  M/M/1  system  with  arrival  rate 
of  Gt  and  to  be  independent  of  each  other. 

The  average  delay  for  link  i  is  given  by  the  waiting  time  in 
the  queue  as 


Where  GC  :.s  the  packet  traffic  in  the  link  and  1/u  the 
average  packet  length  in  bits.  To  form  a  suitable  average 
over  all  the  queuing  processes,  the  T  are  weighed  by  Gc/r 
where  r  is  the  total  input  packet  rate  to  the  network.  In 
this  way,  the  total  delay  suffered  by  all  packets  per  second 
of  network  operation  l  G<- ■  T  is  divided  by  the  total  number 

C 

of  packets  carried  by  the  network  per  second. 

The  total  average  time  delay  per  packet  through  the  net¬ 
work  is  then  given  by 
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Where  Gl/u  =  Ft,  and  L  is  the  total  number  of  links  in 
the  system  and  1/u  is  equal  to  the  average  flow  on  link  . 

It  should  be  noted  that  the  value  of  Cl  (the  minimum  required 
capacity)  must  be  greater  than  the  total  of  all  the  flows  on 
link  i.  In  other  words,  the  value  of  Cl  must  be  greater  than 
the  value  which  is  the  first  run  saturation  level  multiplied 
by  the  arbitrary  capacity  C  for  the  first  run. 

The  value  of  T  includes  both  of  the  waiting  time  and  the 
service  time.  The  waiting  time  is  subject  to  the  interference 
of  all  other  traffic  within  the  network  that  consists  of  the 
data  traffic  as  well  as  the  control  traffic.  This  average 
time  delay  can  be  used  to  analyze  the  relative  goodness  of  the 
routing  strategies  considered  in  this  study. 

In  deriving  this  basic  formula,  a  number  of  factors  (such 
as  processing  time  and  propagation  delay)  are  neglected.  In 
any  realistic  network,  these  variables  as  well  as  others  must 
be  considered  in  the  analysis  of  the  networks.  According  to 
the  above  formula,  the  average  time  delay  depends  only  on  the 
aggregate  flow  of  the  links.  The  reduction  of  the  flow  Ft  of 
each  link  results  in  the  reduction  of  the  average  time  delay 
for  the  network  which  has  constant  input  packet  rate  and  link 
capacities . 

For  the  example  network  (Fig.  3.1),  the  arbitrary  capa¬ 
city  is  10  and  the  value  of  Ct  must  be  greater  than  2.5 
(0.25  x  10).  The  average  time  delay  for  the  capacity  Fl=10, 
and  the  input  (requirement  matrix)  of  2  units  from  node  1  to 
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node  4  and  5  units  from  node  2  to  node  3  are  Tss  =  2.12/r 
for  the  successive  saturation  case  and  Tms  =  2.083/r  for  the 
max-slack  case.  To  make  these  calculations  easy,  the  data 
summarizing  the  final  output  run  are  listed  in  Table  I.  We 
see  that  there  is  not  much  advantage  of  one  procedure  over 
the  other  in  terms  of  performance,  whereas  the  max-slack 
solution  requires  only  two  iterations,  rather  than  6  for  the 
successive  saturation  algorithm. 
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TABLE  I 


i 

! 


T* 

LINK  FLOW  FOR  THE  EXAMPLE  NETWORK 


(  ) ;  Max-Slack 


LINK 

DESTINATION  NODE 

TOTAL 

(Nodes) 

1 

2  3  4 

5 

FLOW 

1-2 

(2) 

(2) 

2 

2 

2-1 

(2.5) 

(2.5) 

2.5 

2.5 

1-3 

(2.5) 

(2.5) 

2.5 

2.5 

3-4 

2-4 

(2.5) 

(2.5) 

2.25 

2.25 

2-5 

0.25 

(2) 

(2) 

2 

2.25 

4-3 

(2.5) 

(2.5) 

2.5 

2.5 

5-4 

0.25 

(2) 

(2) 

2 

2.25 

Req  Mat 

C=10  , 

R]_(4)=2,  R2(3)=5 

T 


III .  EXPERIMENTAL  RESULTS  AND  CONCLUSIONS 

A.  EXPERIMENTS  AND  COMPARISON  OF  THE  RESULTS 

The  experiment  was  conducted  for  two  experimental  networks 
to  analyze  both  approaches,  in  addition  to  the  example  network. 
These  networks  are  specified  in  Appendix  B.  The  networks  were 
used  in  the  experiment  with  requirement  matrices  of  15  units 
from  node  1  destined  for  node  5,  9  units  from  node  7  destined 
for  node  2  for  the  9  nodes/36  links  network  (Fig.  B.l/:  and 
15  units  from  node  2  destined  for  node  12,  4  units  from  node 
5  destined  for  node  6,  and  10  units  from  node  10  destined  for 
node  1  for  the  13  nodes/60  links  network  (Fig.  B.2  .  Both  the 
successive  saturation  and  max-slack  approaches  were  used  for 
the  successive  saturation  approach,  the  21  iterations  are 
required  in  this  experiment. 

The  results  and  their  comparison  is  listed  in  Tables  II, 
III,  and  IV.  These  are  obtained  by  writing  each  flow  variable 
value  from  the  final  output  run  on  each  link  of  the  networks, 
and  summing  these  flow  variable  values  to  obtain  the  total 
flow  on  the  link  (Appendix  F  and  G)  . 

As  shown  in  Tables  II  and  III,  the  number  of  iterations 
required  and  the  total  link  usage  (utilization)  are  much 
greater  for  the  successive  saturation  approach,  but  fewer 
large  link  capacities  are  required  than  in  the  max-slack 
approach.  The  average  time  delay  depends  on  the  capacity 
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TABLE  II 


RESULTS  FOR  THE  13/60  NETWORK 


value  C c  (Table  IV)  .  When  the  capacity  value  CC  is  close 
to  the  minimum  required  capacity,  the  time  delay  is  less  for 
the  successive  saturation  approach  than  for  the  max-slack 
approach  in  both  networks,  and  vice  versa.  It  is  significant 
for  the  large  scale  network,  and  for  large  enough  link  capa¬ 
cities  C t,  that  the  time  delay  difference  becomes  negligible. 

Also,  in  the  process  of  this  study,  it  was  noted  that  both 
approaches  yield  the  same  aggregate  flows  in  the  case  of  very 
simple  networks  (nodes  less  than  or  equal  to  4),  and  for  com¬ 
pletely  symmetric  networks  and  requirement  matrices  (input  and 
capacity  matrices).  For  example,  the  time  delays  are  the  same 
for  both  approaches  with  2  units  input  from  the  outmost  nodes 
to  the  directly  opposite  nodes  respectively  (1  to  13,  2  to  12, 
5  to  9,  12  to  2  and  13  to  1)  for  the  13  node  network,  because 
all  links  total  flow  are  3/4 .  The  detail  experimental  results 
are  listed  in  Appendix  H. 

B.  CONCLUSIONS 

This  is  an  elementary  study  of  the  routing  design  problem 
for  a  packet  network:  given  a  traffic  requirement  matrix, 
minimize  the  average  time  delay  per  packet,  subject  to  finding 
a  feasible  flow  for  a  network  with  fixed  topology  and  link 
capacities . 

The  general  performance  characteristics  and  advantages 
of  the  two  approaches  are  investigated  for  small  and  simple 
networks.  This  illustrates  one  way  to  solve  the  design 
problem  under  limited  conditions. 
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TABLE  IV 


TIME  DELAY  FOR  THE  CAPACITIES 


CAPACITY 

4 

4.2 

4.5  5 

10 

15 

succ 

244. 1/r 

109. 5/r  59. 9/r 

10 . 74/r 

5 . 82/r 

9/36 

SAT 

Net 

MAX 

SLACK 

28. 56/r 

116 . 6/'r  59 . 5/r 

10 . 4 1/r 

5 .7 3/r 

succ 

169. 9/r 

115. 6/r 

75. 7/r  60. 1/r 

15 . 34/r 

3. 9 3/r 

13/60 

SAT 

Net 

MAX 

SLACK 

335. 7/r 

188. 3/r 

114. 2/r  69. 3/r 

14.29/r 

7 . 9 9/r 

APPENDIX  A 

MATRIX  REPRESENTATION  FOR  ILLUSTRATION  NETWORK 
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Figure  A.l.  Matrix  Representation  for  Illustration  Network 
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EXPERIMENTAL  NETWORKS 
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APPEMDI X  C 

MODEL  GENERATION  PROGRAM  FOR  EXAMPLE  PROBLEM 


SJD6  I53SP 

iOPTIONS  FREE 

THIS  IS  THE  DATA  GENERATION  PROGRAM  FuR  MPS^.i. 

THE  FORMAT  ShCJLD  d£  CHANGEu  AOOOROI NGl  Y • 

N  IS  NUMBER  OF  NOCES,  L  IS  NUMBER  OF  LINKS, 
led  »LI  IS  DEFINED  AS  FOLLOwS: 

1.  LINK  INCOMING  TO  NOD  E I  IS  -1. 

2.  LINK  CUTGOING  FROM  NODE  I  IS  l. 

3.  LINK  UNCONNECTED  wITrt  NODE!  IS  0. 

I2CQL<K,D  IS  DEFINED  AS  FOLLOWS: 

1.LINK  OUT GCING  FROM  NODE  K  IS  3. 

2. OTHERWISE,  1. 

I R  (  I  ,  Ki  IS  THE  INPUT  MATRIX  FROM  NODE  I  TO  NODE  K. 
CAP  (  L,  II  IS  THE  LINK  CAPACITY  MATRIX. 

SGLN1! i, UlJ  IS  THE  LINK  MATRICX  CONNECTING  NODES. 


DIMENSION  INC  113  i60l  t  Ic(  15toOI  t  I  ZCDl  <  IS,  oOl  , 

*  IRN(13,6JI  ,  IR(13.13l,CAPlS0'il  .ICOLNil  I  ,o  J  J  , 

*  I C  OLN !  13,601, ICOLN2(I3,oD  I 

INTEGER  I  ,  J , K  ,  L  ? N , N J  L  L , IDEM, K A , J i , I N  C  2 , 1  A , 

*  ID  » Iu , IROnNC, I L ,K0 , KD , IRNI 


*** 
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220 

225 

230 

240 
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21 

23 

20 

19 

0  *** 
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C  **# 


DATA  INITIALl  ZATION  ***** 

READ  !  5»  200  N  ,L 

READ  (5,22  5  )<  I  COLNK1  ,  Jl  )  ,  Jl»l,Ll 
READ (5, 2201  (( IE l  I ,GJ  ,1-1, N1 , J - 1 «  L  J 
READ (5,  2301  U  IRl 1 »K)  , 1= 1 ,NI  , K*1 , Nl 
READ  15,  2  40 )  1C  AP(  I  J,  1  I ,  I  m=  1 ,  L  ) 

FORMAT!  2151 
FORM  AT!  513  J 
FORMAT!  ‘♦110) 

FORMAT!  513) 

FORM  AT! 4E12.5  ) 

ZERO  COLUMN  GENERATION  **** 

00  19  J  =  1,L 

DO  2  0  I  =  l.N 

INCE  =  1 E  1 1 , 0  I 

IF !  I NCE  .EJ.  II  GO  TO  21 

INCE  =  1 

GO  TO  2  3 

INCE  =  0 

IZCQL!  I  »  0  )  =  IN  CE 

CONT INUE 

CONTINUE 

COLUMN  NUMBER  GENERATION  **** 

IAD  =  1000000 
DO  1  K  =  1 ,  N 
I  ROW NO  =  0 
00  12  Ji  =  1 ,  L 

IC0LN!K,J1)  =  ID  0LN1  (  1,  J 1  )  *  K  f  IAD 
IC0L.N2!  K,  Jil  =  ICOLNiK.Jl)  -  K 
CONT INUE 

INCIDENT  MATRICX  GENERATION  *** 

DO  2  I  =  l'N 
IF ( I  .£Q.  Ki  GO  TO  2 
00  3  J  =  1  ,L 
NULl  =  I ZCJLl K, J ) 

IF  ( I  .EG •  KI  NULL  =  0 
IOEN  *  IEU.J  ) 

INC!  I  ,  J  )  *  NULL*  I  OEN 
CONT  INUE 
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IRQ*  NO  =  I « 10  0 

IRMI.R)  *  IRCWNQ  *  K  *  1  AD 

rowT£ecticn  generation  *** 

10=  1 

DO  61  I  -  1»N 
IF  1 1  .EQ.  <)  GO  TO  51 
*RITE<6  ,  140IIRN!  I  »KJ 
CONTINUE 

CAPACITY  ROW  GENERATION  *** 

WRITE  (6  ,  14p)  1C0LN2I1.J) 

coljmnuIection  GENERATION  *** 

WRIT£(6»500i 

00  b5  K  *  1 1 N 

00  o  2  J  =  1  *  L 

NA  =  0 

00  6  3  I  =  i.N 

IF  C X  .EG.  Kl  GO  Tu  53 

NULL  *  I ZCJL (  K»  J  ) 

IDEM  *  I  El  It U  ) 

I NCI  1  i  J  )  *  NU  LL* I DeN 
INC2  =  INCII.J)  ^  Tn 
IFIINC2  .EG.  Qt  GO  TO  6j 

NA  =  NA  ♦  1 

WRITE  (6, 510)  1C0LN  (X,  J  )*  IRNt  I  ,K)  ,  INC-: 

CONTINUE  Tn 

WRITE  (6  *320?  I  COL  N I  Kf  J  ) ,  IC0LN2! K  .  J  >  *  ID 
CJNT INUE 

CONTINUE  „  , 

CAPACITY  GENERATION  *** 

00  66  J  s  1 « L 

IA  =  J  1000  000 

wRI  T  E  ( 6  *  540)  IA*  I  COL  N21 1  *  J)  *  I D 

ALPHA^COLUMN  GENERATION  *** 
wRIT£<6, 565)1  ID 

WRITE  1 6, 570)  ICOLN21  1  ,0  )  t  CAP  I  J ,  1  ) 
CONTINUE 

RHS  GENERATION  *** 

WRITEI6 ,550) 

DO  6  9  K  =  1 1 N 

QQ  6  8  I  *  1  *  N  ^ 

IF 1 1  «E3.  R)  GO  T  0  6  8 
WRIT  £16,560)  IRNU.il, IRCI.RI 
I  CONTINUE 
)  CONTINUE 
00  5  IQ  =  1  ,L 
i  CONTINUE 

WRITE(6,580L  _  _  _ 

)  FORM  ATI  IX,  •  E  C*  ,  17  I 
j  FORMAT!  IX, •  E  L*  *17) 

)  FORMAT!  'COLUMNS'  I 

)  FORM  ATI 4X  X'  •  , 2X,  1  C|  ,  17 , 2X,  I  6 

)  FORMATlix’.'S*  ’,17 ’,1x1  *  L*  J  1 7*  2X»  I  5  ) 

3  FORMAT!  4X*MNPUT  J  ,  5X  ,  JC  •  ,  1 7 , 2X,  I  51 

3  FORM  AT!  •  ENOAT  A*  ) 

STOP 


IC0LN2!  1  ,0  )  .CAP  IJ,  1  ) 


GO  TO  6  8 

IRNI  I,R  ),IR(  I.R  ) 


SENTRY 

5  3 

01  0200 
02  0400 

1-10  0 
1  A  rt 


sentry 


020100 
02050  C 
0 
0 
0 
0 
0 
■1 
0 
1 
0 
0 
0 
0 
0 


-10 

-10 


01  0300 
040300 


030200 

050400 
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APPENDIX  D 


OPTIMIZATION  PROGRAM,  INPUT  jATA 
MATRIX  PICTURE  AND  OUTPUT  FOR  EXAMPLE  PROBLEM 


,/CHANG  JU 

fa&iiPSB*  .  , 

PROGRAM  (  •  NO*  I 

TITL EAV' EXAMPLE  PROBLEM, 
MOVE  t XOB  J, 1  OBJ*  1 

MOVE  (  XRH  St 1  INPUT*  > 

MOVE  (XDATA,'LiNEQS*l 
MOVE  (XPBNAME, 'PROJECT*  > 
CONVERT  ( •  SUMMARY  '  ) 

SETUP  ('MIN') 

3C00UT 

PICTURE 

wrllZARQ 

PRIMAL 

SOLUTION 

EXIT 

P  END 


JOQN  (iCOMU*  ,Ci.A$Ss*a 


SuCC- SAT  *  t 


/* 

//EX5C.SVS  IN 

NAME 

ROUS 

N  OBJ  „  , 

E  C 1 000201 
E  Cl  0003  01 
E  C 1 000401 
E  C 1 00050 1 

e  ci oooi 02 

E  C 1 000302 
E  C1000402 
E  Cl  0005 02 
E  C 1 0001 03 
E  Cl  0002  03 
E  C100Q403 
E  C1000503 
E  Cl  0001 04 
E  Cl  0002 04 
E  C1000304 
c  Cl  0005 04 
E  Cl  0001  05 
E  Cl  0002 05 
E  Cl  0003 05 
6  Cl  000405 
c  Ll  0102  00 
E  Ll 0201 00 
E  Ll 010300 
E  Ll  030200 
E  L 1020400 
E  Li  0205  00 
E  L 10403  00 
E  L 1050400 
COLUMNS 

X1020i  01 
XlQ2Qi 0 1 
X 1 0302  01 
X1030201 
XI 0302  01 
XI 02040 1 
X 1 02040 1 
XI 020401 


00  * 
LINEQS 


C 1000  201 
L 1020  100 
C 1 000  201 
C 1 000 301 
L 1030  200 
Cl 000  201 
C 1000  401 
Ll 020  400 


I 

1 

-1 

1 

1 

1 

-1 

1 


XI 020501 

C 1000  201 

1 

XI 0205  01 

C 1 000  501 

-1 

Xl 0205  01 

11020500 

1 

X1040301 

C1000301 

-1 

XI 0403 01 

C1000401 

1 

XI 0403 01 

L 1040 200 

1 

XI 050401 

C1000401 

-1 

Xl 050401 

C1000501 

1 

Xl 050401 

L 1050  400 

1 

Xl 0102  02 

C1000  102 

1 

X101 02  02 

L 1 01 0200 

1 

Xl 010302 

C1000102 

1 

Xl 010302 

Cl  000  302 

-1 

Xl 010302 

L 1010200 

1 

X1030202 

Cl  000  302 

1 

X 1 0302  02 

L1030200 

1 

Xl 040302 

C 1000  202 

-1 

X 10403  02 

C1000402 

1 

Xl 0403 02 

L1040200 

1 

X1050402 

Cl 000 402 

-1 

Xl  050402 

C1000  502 

1 

Xl 050402 

L 1050 400 

1 

Xl  010205 

C1000  103 

1 

Xl 0102  03 

C1000203 

-1 

Xl 0102  03 

L 1010  200 

1 

Xl 0201 03 

C1000  103 

-1 

X1020103 

C 1000203 

1 

Xl 0201 03 

L1020100 

1 

X 1 0103  03 

Cl  000103 

1 

X 1 0103  03 

L 10 10  200 

1 

X 1 020403 

Cl 000  203 

1 

Xl 020403 

C1000403 

-1 

X 1 020403 

L102Q 400 

1 

Xl 0205 03 

C 1000203 

1 

Xl 0205  03 

C 1000  503 

-1 

X 1 0205  03 

L1020  500 

1 

Xl 0403 03 

C 1000  403 

1 

X 10403 03 

L 1040 200 

1 

X 1 050403 

Cl  000 403 

-1 

X 1 050403 

Cl  000  503 

1 

Xl 050403 

(.1050  400 

1 

Xl  0102  04 

C 1000 104 

1 

Xl 0102  04 

Cl 000 204 

-1 

Xl 0102  04 

L10102Q0 

1 

X 10201  04 

Cl  000  104 

-1 

X 1 020104 

C 1000  204 

1 

Xl 0201 04 

Ll  020  100 

1 

X1010304 

C1000104 

1 

Xl 0103  04 

C1000  204 

-1 

Xl 0103  04 

Ll 010200 

1 

Xl 0302 04 

C 1000 204 

-1 

Xl 0302  04 

C1J00304 

1 

Xl 0302  04 

Ll 030  200 

1 

X 1 020404 

C 1000 204 

1 

Xl 020404 

Ll  02  0  400 

1 

Xl 020504 

C 1000  204 

1 

X 1 020504 

C 1000 £04 

-1 

Xl  0205  04 

L1020  500 

1 

Xl  05  04  04 

Cl 000  £04 

1 

Xl 050404 

L1050400 

1 

X 1 0102  05 

Cl  000  105 

1 

Xl 010205 

CiOOO  205 

-1 

X 1 01 02  05 

Ll 010  200 

1 

X 1 0201 05 

C1000105 

-1 

X 1 0201 05 

CIOOO  205 

1 

X 1 0201 05 

L102Q100 

1 

Xl 010305 

CIOOO  105 

1 
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xioiosos 

C1J00305 

-1 

Xl  0103  05 

L 1010  300 

1 

X 1 0302  05 

C10Q0205 

_ > 

Xl 0302  05 

C1000305 

i 

X 1 030205 

L 1030 200 

1 

X 1 020405 

C 1000  205 

1 

X 1 02040  5 

C1000405 

-1 

X 1 020405 

L1020400 

1 

Xl 0205  05 

C 1000  205 

1 

Xl 0205  05 

L 1020 500 

1 

Xl 0403 05 

Cl 000 305 

-1 

Xl 040305 

C 1000 405 

1 

Xl 0403  05 

L 1040 300 

1 

S1000001 

L1010200 

1 

SI  0000  02 

L 1020100 

1 

SI  000003 

Ll 010300 

1 

SI  0000  04 

L 1030  200 

1 

SI  000005 

L 1020  400 

1 

SI  000006 

Ll 020 500 

1 

SI  0000  0  7 

Ll 040  300 

1 

SI  0000 08 

1 1050  400 

1 

AL  PHA 

OBJ 

1 

ALPHA 

L 1010  200 

-0. 10000E 

02 

ALPHA 

L1020100 

-0.10000E 

02 

ALPHA 

Ll 010  300 

-0.10000c 

02 

ALPHA 

L 1030  200 

-0. 10000E 

02 

ALPHA 

Ll 020  400 

-0 • 10000c 

02 

ALPHA 

Ll 020  500 

-O.iOOOOE 

02 

ALPHA 

L 1040 300 

-0.10000E 

02 

ALPHA 

L1050400 

-0 .10  000E 

02 

INPUT 

C1000  201 

0 

INPUT 

C1000  301 

0 

INPUT 

C 1000 4Q1 

0 

INPUT 

C1000501 

0 

INPUT 

Cl  000  102 

0 

INPUT 

C 1000  302 

0 

INPUT 

C1000402 

0 

INPUT 

C1000S02 

0 

INPUT 

C1000103 

0 

INPUT 

C1000  203 

5 

INPUT 

C 1000403 

0 

INPUT 

C1000503 

0 

INPUT 

Cl 000  104 

2 

INPUT 

C 1000  204 

0 

INPUT 

Cl  000  304 

0 

INPUT 

Cl  000  504 

0 

INPUT 

C1000 105 

0 

INPUT 

C1000  205 

0 

INPUT 

C 1000  305 

0 

INPUT 

C1000405 

0 

ENQATA 

/* 

// 
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APPENDIX  E 


OPTIMIZATION  PROGRAM,  INPUT  DATA 
MATRIX  PICTURE  AND  OUTPUT  FOR  MAX-SLACK  APPROACH 


//CHANG  JO  8  < 1538,  1808),  'CHANG  JOQN  rfOGNS •  ,CLASS  =  A 

//* 

//*  MPSIII  EXAMPLE  PROBLEM 

//* 

//  EXEC  MSSMPS 
//CPC.SYSIN  DO  * 

PROGRAM  ( •  ND*  ) 

I NITIALZ 

TITLE  (‘EXAMPLE  PROBLEM,  END  RUN  FOR  MAX-SLACK.') 
MOVE  (XOB J, ' OBJ'  ) 

MOVE  (  XRH  S, '  I  NPJ T * ) 

MOVE  ( XDA  TA»  •  LI  \l  £  OS  '  ) 

MOVE  (XPSNAME,' PROJECT*  ) 

CONVERT  (  '  SUMMARY  •  ) 

SETJP  ('MAX') 

3  COQUT 

PICTURE 

KriUA  RC 

PRIMAL 

SOLJTION 

EXIT 

PEND 


/* 

//EXEC.SYSIN 
NAME 
ROmS 
N  OBJ 
E  C1000201 
E  C 1 0003  01 
E  Cl  000401 
E  C 1 0005  01 
E  C 10001 02 
E  C 1 0003  02 
E  C 1 000402 
E  C 1 0005  02 
Cl  0001 03 
Cl  0002 33 
Cl  000403 
Cl  0005 03 
C 1000104 
Cl  0002 04 
C 1 000304 
C 1 000504 
C 1 0001 05 
C 1 0002  05 
Cl  0003  05 
C 1 000405 
L1010200 
LI  0201 00 
L1010300 
L 1030200 
LI  0204 00 
L 1 0205  00 
LI  0403 00 
L 1 0504  00 
CQLJMNS 

XI 0201 01 
X 1 0201 01 
X 1 030201 
Xl 0302  01 
Xl 030201 
X 1 020401 
Xl 020401 


c 

E 

E 

E 

E 

E 

E 

1 

E 

E 

c 


E 

E 

E 

E 


CD  * 
LINEQS 


C1000201 
LI  020 100 
C 1000  201 
C1000  301 
L1030200 
C1000  201 
C 1000  401 


1 

1 

-1 

1 

1 

1 

-1 


u 


X 1 020401 

Li 020 400 

1 

Xl 0205  01 

C 1000  201 

1 

Xl 020501 

C 1000  501 

-1 

Xl 0205  01 

Ll 020  500 

1 

Xl 0403  01 

C1000  201 

-1 

X 10*03 01 

C 1000 401 

1 

Xl 0*03 01 

L 1040 300 

1 

Xl 050401 

Cl  000  401 

-1 

Xl 050401 

Cl  000  501 

1 

X1050401 

L 1050 400 

1 

Xl 010202 

C 1000 102 

1 

Xl 0102  02 

L 1010 200 

1 

Xl 010302 

C 1000102 

1 

Xl 010302 

C 1000 502 

-1 

Xl  010302 

L1010300 

1 

X1030202 

C 1000  202 

l 

Xl 0302  02 

L1030200 

1 

Xl  0403  02 

Cl 000  302 

-1 

X 10403 02 

C 1000  402 

1 

Xl 040302 

L1040 300 

1 

Xl 050402 

C 1000402 

-i 

Xl 050402 

C 1000 502 

1 

X 1 050402 

L1050400 

1 

X1010203 

C1000103 

1 

Xl 0102  03 

C1000  203 

-1 

Xl  010203 

L1010200 

1 

X 10201 03 

C1000  103 

-1 

Xl 0201 03 

C1Q0Q202 

1 

X1020103 

L 1020 100 

1 

Xl 0103  03 

C 1000 103 

1 

Xl 0103  03 

L1010300 

1 

X 1020403 

Cl  000  203 

1 

Xl 020403 

Cl  000 403 

-1 

X 1 020403 

L 1020  400 

1 

Xl 0205  03 

C 1000 203 

1 

Xl 0205  03 

Cl  000  503 

-1 

Xl 020503 

Ll 020  500 

1 

Xl 0403  03 

C  1000  403 

i 

X 1 0403  03 

L 1040 200 

1 

X 1 050403 

C 1000403 

-1 

Xl 050403 

C1000  503 

1 

Xl 050403 

L 1050  400 

1 

Xl 0102  04 

C1000  104 

1 

X 1 0102  04 

C1000 204 

-1 

X 10102 04 

L 101 0200 

1 

X 1 0201 04 

C 1000  104 

-1 

Xl 0201 04 

C1000204 

1 

X 1  0201 04 

L1020100 

1 

Xl 0103  04 

C1000  104 

1 

Xl  010304 

Cl  000 304 

-1 

Xl 010304 

L 1010  200 

1 

Xl 030204 

Cl 000 204 

-1 

X1030204 

Cl  000 204 

1 

X 1 030204 

L 1030 200 

1 

Xl 020404 

C 1000  204 

1 

Xl 020404 

L 1020  400 

l 

X 1 0205  04 

C1000  204 

1 

Xl 020504 

Cl  000 504 

-1 

Xl 020504 

L 1020  500 

1 

X 1  0504  04 

CIOOO 504 

1 

Xl  0504  04 

L1050400 

1 

X101Q205 

C1000  105 

1 

Xl 0102  05 

Cl  000 205 

-1 

X 1  0102  05 

L 1010 200 

1 

Xl 0201 0  5 

C1000 105 

-1 

Xl 0201 05 

C 1 000  205 

1 

Xl 0201 05 

L 1020  100 

1 
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Tm 


r 

[ 

XI 0103  05 

C1J00105 

1 

. 1 

! 

XI 0103  05 

C 1000 205 

-1 

XI 0103  05 

21010300 

1 

XI 0302 05 

C 1 000  205 

-1 

\ 

X 1 0302  0  5 

Cl  000 305 

1 

XI  030205 

L 1030  200 

1 

! 

X 1 020405 

C 1000  205 

1 

XI 020405 

C1000405 

-1 

X 1 020405 

L1J20400 

1 

XI 0205  05 

C 1000  205 

1 

XI 0205 05 

L1020500 

1 

XI 0403  05 

C1000205 

_  1 

XI 0403  05 

Cl  000 405 

I 

XI 0403 05 

LI  040300 

1 

SI  0000 01 

060 

1 

t-i  01  020  0 

1 

Si  000002 

33  J 

1 

Ll 020i0  0 

l 

SI  0000 03 

OBJ 

1 

L1010300 

1 

S 1 0000  04 

OBJ 

1 

L1030200 

l 

SI  000005 

OBJ 

1 

L 102040  0 

l 

S 1 0000  06 

OBJ 

1 

l1Q2050  0 

l 

SI  000007 

OBJ 

1 

Li  04030  0 

l 

SI  0000  08 

OBJ 

1 

Li  05040  0 

i 

t 

AL  PH  A 

OBJ 

1 

* 

ALPHA 

L1010200 

-0. 10000E 

02 

♦ 

ALPHA 

L 1020  100 

-0  •  IQGOQE 

02 

* 

ALPHA 

LI  01 0300 

-0 .100006 

02 

* 

ALPHA 

L1030200 

-0. lOOOOE 

02 

* 

ALPHA 

L1020400 

-0.10000E 

02 

* 

ALPHA 

LI  020  500 

-0 .lOOOOE 

02 

* 

ALPHA 

LI 040 300 

-0. lOOOOE 

02 

* 

AL  PHA 

LI  050  400 

-0. lOOOOE 

02 

*HS 

INPUT 

C 1000  201 

0 

INPUT 

C1000  301 

0 

INPUT 

C1000401 

0 

INPUT 

C 1000  501 

0 

INPUT 

C1J00  102 

0 

INPUT 

C 1000  202 

0 

INPUT 

C1000  402 

0 

INPUT 

C 1 000  502 

0 

INPUT 

C1000  103 

0 

INPUT 

C1000  203 

5 

INPUT 

Cl  000 403 

0 

INPUT 

C  1000  503 

0 

1 

INPUT 

C100Q104 

2 

INPUT 

C 1 000  204 

0 

INPUT 

C1000304 

0 

INPUT 

Cl  000 504 

0 

INPUT 

C1000  105 

0 

INPUT 

C 1 000  205 

0 

INPUT 

C1000205 

0 

I N  PUT 

C 1000 405 

0 

INPUT 

L1010200 

2.5 

INPUT 

L1020  100 

2.5 

INPUT 

LI  01 0  300 

2.5 

INPUT 

Ll  030  200 

2.5 

! ' 

INPUT 

L 1020 400 

2.5 

INPUT 

L1020500 

2.5 

|  ; 

INPUT 

Li  040 200 

2.5 

i 

INPUT 

L 1050 400 

2.5 

*  ■ 

1  ENOATA 

i : 

i  >* 

,  j 

\  // 

i 

i  * 
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HUJi 


DESTINATION  NODE 


TOTAL 


(NODES)  1  2  3 


REQ-MAT 


^ (5)=15,  R5(2)=9 
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LINK  FLOW  FOR  13/60  NETWORK 


9-2 


LINK  DESTINATION  NODE 


_ DESTINATION  NODE _ TOTAL 

(NODES)  1  2  3  4  5  6  7  8  9  10  11  12  13  FLOW 


LINK  DESTINATION  NODE  TOTAL 


LINK  FLOW  FOR  SYMMETRIC  REQ-MAT  OF  13/60  NETWORK 
Table  H.l.  Link  Flow  for  Symmetric  Req-Mat  of  13/60  Network  (  );MAX-SLACK 


DESTINATION  NODE 


DESTINATION  NODE 


LINK  OESTINAl ION  NODE 


LIST  OF  REFERENCES 


1.  Arne,  T. ,  Computer  Methods  in  Operations  Research,  p.  79, 
Academic  Press ,  1978 . 

2.  Pooch,  U.  W. ,  Green,  W.  H.,  and  Moss,  G.  G. ,  Telecommuni¬ 
cations  and  Networking,  Little,  Brown  Computer  Systems 
Series,  1982. 

3.  Gass,  S.  I.,  Linear  Programming  Method  and  Application, 

4th  Ed.,  McGraw-Hill,  1975. 

4.  Kleinrock,  L. ,  Queueing  Systems,  V.  2,  pp.  292-297,  Wiley, 

19  76  . 

5.  Davices,  D.  W. ,  and  Barber,  D.  A.,  Communication  Networks 
for  Computer,  pp.  291-294,  Wiley,  1973. 

6.  Carvis,  H.,  Communications  Network  Analysis,  pp.  103-104, 
Lexington,  1981. 

7.  Hiroshi,  I.  and  Tadao,  S.,  Theoretical  Aspects  in  the 
Analysis  and  Synthesis  of  packet  Communication  Networks, 
Proceedings  of  the  IEEE ,  V.  66 ,  No .  11,  No vT§75"I 

3.  Kuo,  F.  F.,  Protocols  and  Techniques  for  Data  Communications 
Networks ,  pp.  124-127,  Prentice-Hall ,  1981. 

9.  Ketron  Inc.,  MPS  III  User  Manual,  1980. 

10.  Phillips,  D.T.,  and  Alberto,  G.D.,  Fundamentals  of  Network 
Analysis ,  Prentice  Hall,  1981. 


33 


INITIAL  DISTRIBUTION  LIST 


No . 


Defense  Technical  Information  Center 
Cameron  Station 
Alexandria,  VA  22314 

Library,  Code  0142 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  John  W.  Wozencraft,  Code  6  2VJz 
Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Professor  Paul  H.  Moose,  Code  62Me 
Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Lieutenant  Colonel  Chang  Joon  Wong 
Korea  MC 

Navy  Headquarters  (N-l) 

Seoul  Korea 


Copies 

2 

2 

3 

2 

3 


84 


